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Abstract 

In this paper wc consider a jump-diffusion dynamic whose parameters are driven by a 
continuous time and stationary Markov Chain on a finite state space as a model for the 
underlying of European contingent claims. For this class of processes wc firstly outline the 
Fourier transform method both in log-price and log-strike to efficiently calculate the value of 
various types of options and as a concrete example of application, we present some numerical 
results within a two-state regime switching version of the Merton jump-diffusion model. Then 
we develop a closed-form solution to the problem of pricing a Forward Starting Option and 
use this result to approximate the value of such a derivative in a general stochastic volatility 
framework. 

Key words: regime switching jump-diffusion models, option pricing, Fourier transform 
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1 Introduction 

Since the paper by Naik (1993), the use of continuous time regime-switching processes to model 
asset price dynamics stimulated an increasing interest in the context of option pricing. The 
empirical evidence of a regime switching behavior of some economic time series was pointed 
out by Hamilton (1989, 1990), who suggested the use of an underlying Markov chain switching 
between regimes to account for some peculiarities in observed data. The ability of these econo- 
metric models to capture specihc features such as volatility clustering and structural breaks is 
widely recognized (see e.g. Timmermann (2000)). Consequently, they can be considered as an 
appealing class of models also in the framework of derivative pricing. In the last decades there 
has been a considerable progress in the pricing exercise for plain vanilla European or American 
style options: see e.g. Di Masi et al. (1994), Bollen (1998), Guo (2001), Hardy (2001), Duan et 
al. (2002), Buffington and Elliott (2002), Konikov and Madan (2002), Guo and Zhang (2004), 
Chourdakis (2004,2007), Edwards (2005), Liu et al. (2006), Yao et al. (2006), Jobert and Rogers 
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(2006), Elliott and Osakwe (2006), Jiang and Pistorius (2008), Boyarchenko and Levendorskii 
(2009), Khaliq and Liu (2009), Di Graziano and Rogers (2009), Ramponi (2009), Liu (2010). 
Comparatively few results are available for exotic options: see Boyle and Draviam (2007) and 
Elliott et al. (2007). Such results typically differ in the model considered (switching diffusions or 
more general Levy processes), in the technique for solving the pricing problem (direct evaluation 
of expectations with respect to the probability density of the underlying, numerical solution of 
the associated PDE, recombining trees, Fourier transform methods) or in the type of financial 
product to price. 

In this paper we consider a quite general underlying dynamic which can be seen as a switching 
Levy process of the I type, or finite activity Levy process (see e.g. Cont and Tankov (2004)). 
In particular, on a filtered probability space (fi, {J^t}, J-", P) the dynamic is of the form 

S{t) = soe^(*) (1) 

where X{t) is specified as a jump-diffusion whose parameters are driven by a{t), a continuous 
time and stationary Markov Chain on the state space 5 = {1, 2, . . . , M}. This model provides an 
example of non-affine and non-Levy process for which we are able to calculate the characteristic 
function (see Prop. 13. ip and therefore the pricing problem for European style options is efficiently 
faced through Fourier transform techniques. Such techniques, originated by the works of Heston 
(1993) and Carr and Madan (1999), are based on the representation of the value of the option 
in a proper Fourier space and have been successfully applied to a variety of pricing problems in 
the last years. Among the various contributions to this theory, see Bakshi and Madan (2000), 
Raible (2000), Lewis (2002), Lee (2004), Hubalek et al. (2006), Biagini et al. (2008), and more 
recently Cherubini et al. (2009), Dufresne et al. (2009), Hurd and Zhou (2010), Eberlein et 
al. (2010). Following this approach we can price various types of European options under the 
regime-switching dynamic by using the Fourier transform method both in the log-price space and 
in the log-strike space, consequently taking advantages from the powerful Fast Fourier Transform 
(EFT) computational tool. The case for a switching pure jump process has been considered in 
Ehiott and Osakwe (2006). 

As an application we consider the problem of valuing a Forward Starting option (FSO) for 
which an almost (i.e. up to numerical integration) closed-form solution is obtained in term of 
an integral transform. A similar technique was used in Kruse and Nogel (2005) to price a FSO 
in the Heston stochastic volatility model. These options are well-known exotic derivatives (see 
e.g. Hull (2009)) characterized by the payoff 

UT{S{T),K) = S{T)-KS{t*) (2) 

where t* S (0,T) is the determination time and k € (0, 1) is a given percentage. They are the 
building blocks of the so-called cliquet options. As it will be shown, our formula is very simple, 
being a finite mixture of call prices evaluated at the determination time under each regime, 
weighted by the stationary probability of the chain. Furthermore, in Chourdakis (2004) a pro- 
cedure to approximate the value of an European option in a model with stochastic volatility and 
jumps was proposed by building a continuous-time Markov chain which "mimics" the volatility 
process. The approximating dynamics turns out to be a regime-switching jump-diffusion model. 
By using such an approximation, a pricing algorithm for FSO in a general stochastic volatility 
model can be designed based on our mixture representation. 
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The paper is organized as follows. In Section 2 the dynamic model for the underlying is 
presented with a scheme for its numerical simulation and a useful representation through the 
sojourn times of the underlying Markov chain is introduced. In Section 3 the Fourier transform 
method both in log-price and log-strike is considered and formulas for the price of an European 
call option are explicitly derived. A numerical example of calibration on real data for a two-state 
regime switching jump diffusion model with gaussian jumps is reported. Finally, in Section 4 
the price of a Forward starting option is obtained by using the Fourier transform representation 
and an algorithm to get approximate prices in a general stochastic volatility model is outlined. 

2 The model 

Let us consider on a filtered probability space {0,,{Tt},J-,V), the asset price dynamic of the 
form 

Sit) = soe^(*) (3) 

where X{t) is specified as follows. 

Let a(t) be a continuous time, homogeneous and stationary Markov Chain on the state 
space 5 = {1, 2, ... , M} with a generator Q G M*"^^*'^; furthemore ^ : S ^ R, a : S ^ R and 
J : E X S ^ M. are given functions, {E, £) being a measurable mark space. In a given interval 
< t < T, we consider the following dynamic 

dX{t) = i{a{t))dt + a{a{t))dW{t) + dJ{t), X(0) = 0, 

J{t) = f I j{y,a{s-))p''{dy,ds) 
Jo Je 

where p" {dy , ds) is a marked point process (Runggaldier (2003)) characterized by the intensity 

X{a,dy) = \{a)m{a,dy). 

Here A(-) represents the (regime-switching) intensity of the Poisson process Nt{E), while m(-, dy) 
are a set of probability measures on {E,£), one for each state (regime) « G 5 of the chain. The 
function 7(2/, a) represents the jump amplitude relative to the mark y in regime a. Throughout 
the paper we assume that the processes a(-) and W{-) are independent and that W{-) and 
p°'{dy,dt) are conditionally independent given a{t). We denote = a{a{s) : < s < t} 
the (T-algebra generated by the Markov chain. Furthermore, we assume that E[eT(^(°)'°)] = 
e"'^^'"^ m{a, dy) is finite for each regime a, where Y{a) is the random variable associated to 
the measure m{a,dy). We also define the compensated point process q°{dy,dt) = p°'{dy,dt) — 
X{a{t—))m{a{t—),dy)dt in such a way 

f f H{y,a{s-))q"{dy,ds) 
Jo Je 

is a martingale in t for each predictable process H satisfying appropriate integrability conditions. 
In particular the jump process 

Jit) = l{y,a{s-))q°'{dy,ds) + I / -i{y,a{s))\{a{s))m{a{s),dy)ds 

Jo Je Jo Je 

= J{t)+ f A(a(s))E[7(y(a(s)),a(s))]ds, Y {a{s)) ^ m{a{s) , dy) 
Jo 
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is the sum of a martingale and an absolutely continuous process, whenever 7 satisfies the proper 
conditions. 

A sample path of this process is generated as follows (see Fig.l): 

1. generate a path of the Markov chain, i.e. a set of switching times tq = 0, ti, . . . , r^, r^^i = 
T and the corresponding states a{t) = at ^ S,Tk <t < Tfc+i, A; = 0, . . . , L; 

2. generate the jump times Vk^ of the Poisson process in each interval [rfc,r/c+i) according to 
the intensity X{ak) and let Nk be the number of jumps; 

3. for any k = 0,...,L generate Nk i.i.d. samples Y{a{vkj)),j = l,...,Nk distributed 
according to the probability m{ak,dy); 

4. on a given time grid to, . . . ,tn of [0, T] built as the superposition of a deterministic grid 
and the jump times Vj, let X(to) = and 

= X{t,)+^{a{ti)){t,+^-ti)+a{a{ti)){W{ti+i-W{t,)), (4) 

X{ti+i) = X{ti+i-)+ -f{y,a{ti+i))p{dy,ti+i). (5) 

Je 

If tj+i is actually a point of the Poisson random measure, the magnitude of the jump is 
sampled, that is 

/ liy,aiti+i)p{dy,ti+i) = j{Y{a{ti+i)),a{ti+i)), 
Je 

otherwise the jump term is zero. 



In view of our pricing application, from now on we assume to specify our model in a prob- 
ability space where the "discounted" asset price is a martingale. In 
particular, we keep the function ^ : 5 — )■ M unspecified in order to cope with slightly different 
types of contracts: for example, given the risk-free rate r, we can set n{a) = fi with fi = r — q, 
q being the (continuous) dividend rate, ^ = r — rj, rj being the foreign risk- free rate or more 
generally ii{ot) = r{a) — q{a) if rates and dividend are regime-switching too. 

We therefore consider the following model 



1 



X{t) = j l^fi{a{s)) - ^a^{a{s)) - X{a{s))K{a{s))j ds + J a{a{s))dW{s) 

+ f I 7{y,a{s-))p-{dy,ds), (6) 



E 



where 



K(a) =E[(e^(^(")'") - 1)], a£S. 
An application of the generalized Ito's Formula gives 



dS{t) = S{t-) i^-X{a{t))K{a{t))dt + a{a{t))dW{t) + J (e^(?^'"(*-» - l)p'^{dy,dt) 
= S{t-) (a{a{t))dW{t) + / (eT(s^'°(*-)) - l)q''{dy,dt)] , 



E 
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Figure 1: A sample path of the RSJD - Example l2.1l Red cross are the switching times, circles represent 
the jump times for the two regimes. 



where q°'{dy,dt) = p"{dy,dt) — X{a{t—))K{a{t—))dt is the compensated process. Hence, S{t) is 
a martingale. The corresponding jump-diffusion SDE for the asset price is therefore 



= {fi{a{t))-\ia{t))K{amdt + a{a{t))dW{t) + |^(e^(^'°(*-)) - dt), 5(0) 



So- 

(7) 



Example 2.1 As a working example we consider a two-state regime switching version of the 
Merton jump-diffusion model. This is defined by taking 'y{y, a) = y and two kinds of normal 
jumps, i.e. Y{i) ~ J\f{ai,bi) from which K{i) = E[(e^(*) — 1)] = e'^'+''»/^ — 1, i = 1,2. The two 

state Markov chain a{t) € 5 = {1,2} has generator Q = [ ). Let ai,Xi > and 

V Q2 -Q2 ) 

fii,i = 1,2 be given parameters: the regime switching jump- diffusion Merton model is defined as 
dX{t) = [n{a{t)) - ^a^{a{t) - X{a{t)) K{a{t))]dt + a{a{t))dW{t) + dJ{t) 



Jit) = / / yp'^idy,ds), X{t,a{t),dy) = X{a{t))(pa{t)iy)dy 
Jo Je 

where X{a{t)) G {Ai,A2}, a{a{t)) G {cti,cj2}, fi{a{t)) G {/xi,/i2} and X{t,a{t),dy) is the intensity 
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process of the Poisson jump component, (j)i{y) being the density of a normal distribution N {ai, bi), 
z = 1,2. □ 

Next Proposition gives a useful representation for X{T). A sketch of the proof is reported 
in the Appendix. 



Proposition 2.1 Let Ti = I^^^^^jds, i = 1,...,M be the occupation times of the Markov 

rT A/ 



chain and let us define ^(a) = ir{a) — icr^(a) — \{a)K{a) and 



■t{Ti,...,Tm)= [ ^{a{s))ds = f2^i^Ti. 

■^0 i=l 



Then the process X(T) admits the following representation: 

M M N{Ai) 

X(r) = HHTi,...,r^) + ^a(i)Z(A,) + ^ Y^, (8) 



i=l i=l k=l 



where N{Ai) and Z{Ai) are distributed as Poisson variables Poiss(AjTj) and as Normal variables 
Af{0, Ti), respectively, i = 1, . . . , M . 



It is readly seen that by defining Xt^x = X{T) — X{t) we have 

M M N{A,) 

X,T = H,T(T*'^,...,Tl,f)+J]aWZ(A,)+J] (9) 

1=1 1=1 k=l 

where T^''^ = I^(^g^^jds, j = 1,...,M, and the random variables N{Ai) and Z{Ai) have 
conditional distributions Poiss(Aj(rj — t)) and Af{0, Ti —t). Correspondingly we can write S(T) = 
S{t)eMXt,T). 

Remark 2.1 Notice that for a Levy process XL{t) having characteristic function 

E[e-^i(*)] = exp {t(iiu - ^ + ^(e'"" - lWx))j 
the expected value is 'Ei[Xi{T)] = T[^ + xf3{dx)). For our RS model, it follows from \2. 1\) that 

E[X(r)] = Yi^ii) + A(i)/Zy(o)E[ri] = Yi^ii) + A(i)^y„) / V{a{s) = i\a{0))ds. 
i=i i=i 

This quantity can be easily evaluated for a two-state MC, since we have Ps = e*^^ = 
^ ( 1.(1 _ e-^(M+^)) ue-'^^'+''^ + fi )' ™^ ^^nplies that, starting e.g. from a{0) = 1 

mm] = +A(iw..,)(^ +^r) 

+ (e(2) + A(2K(.))(-^ ^— + ^T). 

jJL + U IJL-\- V 11 + V 
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3 The transform method 

Our main interest is the efficient numerical evaluation of the price Ho of an European contingent 
claim specified by the payoff function n(s, K), exercised at the future time T, K being a trigger 
parameter. By letting r(t) be the interest rate process, B{t) = eiip{fQ r{u)du) the usual money 
market account and P{t, T) the time-t value of a discount bond maturing at T, arbitrage pricing 
theory and a change-of-numeraire technique give the well-known characterization of prices 

Ho = B'^[B{T)-^U{SiT),K)] = P{0,T)'E^[U{S{T), K)], P{0,T) = E^[B{T)-^], 

where V is the risk-neutral measure and Q is known as T- forward measure. When interest rates 
are deterministic, the two measures are equal. In the following we assume that our dynamic 
model is given under the measure Q. All the expected values will be considered with respect to 
this measure. 



It is well known that Fourier transform methods can be efficiently used for the valuation of 
European style options. Two main variants have been developed depending on which variable 
of the payoff is transformed into the Fourier space. In view of the structure assumed for the 
dynamic of the underlying price S{T) and our next applications, we consider instead log(S'(T)) = 
X{T) + log(so) as the state variable and k = log(i^) for the trigger parameter, in such a way 
for any payoff Il{s,K) = n(e^°^'^^\ 6^°^^^^-') = Il{y,k). Correspondingly, we can consider the 
generalized Fourier transform with respect to the state variable y, Ilfc(z) = f^e^^'^Il{y, k)dy 
(log-price transform), or w.r.t. the trigger k, 'ily{z) = f^e'^^''Il{y,k)dk (log-strike transform), 
z € C. In general we assume that these transforms exist in some strip Sji = {2; € C : —00 < 
a < Q{z) < b < -|-oo}0 of the complex plane. Examples of payoffs are reported in Table 
The first approach was proposed in this form in Raible (2000) (but the representation of option 
prices through inversion of characteristic function appeared for the first time in Heston (1993)), 
while the second was introduced in Carr and Madan (1999). 

Formally, Fourier inversion gives 



U{y,k) 




e ^^yilk{z)dz, log-price transform 



e Ily{z)dz, log-strike transform 



where integrals are considered along the straight line '^{z) = v in the complex plane. By letting 
^Tiz) = E[e^^''''''"^-*], z G C be the (generalized) Fourier transform (or characteristic function) of 
X{T), we have 

no/P(0,r) = E2[n(X(r)+log(so),A:)] = / U{y,k)QT{dy) 



Li^i::-::^-'''mz)dzQT{dy) 



1 f iv+00 ■ 

2tt 



£-::^k{z)f^e-yQT{dy)d2 



Q{z) and ^{z) stand for the imaginary and real part of a complex number, 2 — ^{z) + i^{z) £ C. 
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Payoff 


GFT in log-price 


Strip of regularity 


GFT in log-strike 


Strip of regularity 


{ey - e^)+ 


gfcCiz + l) 


> 1 


iz—z'^ 


< 


(e*^ - ey)+ 


gfe(iz + l) 

iz—z'^ 


< 


e!/(iz+l) 
iz— 2"^ 


> 1 




a+iz 


> a 


12 


> 


min(e^, e'^) 


z'-'—iz 


< < 1 


z^—iz 


< < 1 



Table 1: Generalized Fourier transforms of typical payoffs. 



'i^ liu-x! ^ (— log-price transform 

^/ir-(^^~'''''^^[nx(T)+iog(so)(^)]^^' log-strike transform. 

In order to justify the previous equalities, some conditions are required: existence of the 
generalized Fourier transform 11, integrability along the contour Si(-z) = in some strip 5n in 
order to guarantee the Inversion Theorem and existence of the expectation E®[e''^("^^] (see Lee 
(2004) for log-strike transform, Lewis (2002) or the recent Eberlein et al. (2009) for log-price 
transform). Notice that the use of generalized Fourier transform permits to exploit contour 
variations by means of the residue theorem, as it will be seen in next paragraphs. 

Due to the exponential structure of the GFT of typical payoffs (see Table ([I])), also for the 
log-strike transform it is required the calculation of '■Pt{z) appearing through the expectation 
'E,'^\Ylx{T)+\og(so){^)]- Next Proposition gives the GFT of our process. Similar results are avail- 
able (see Chourdakis (2004)) where a particular structure of the generator Q is considered: for 
completeness, we report the proof in the Appendix. 

Proposition 3.1 Let (j)j{z) = 'E[e^^^^'^^^^'^^] be the generalized Fourier transform of the jump 
magnitude. Then, by letting 

M^) = ^e(j) + ^i^V2(i) - iX{j)iMz) - 1) (11) 

and ^i{z) = ^j{z) — i^Miz), we have 

iprp(^z) = e''^M{z)T (^1 . diag(^i(2),...,,9M-l(2),0))T . J(^q^^ 



(12) 



where 1 = (1, . . . , 1)' G M^^^^, 1(0) = (I„(o)=i, • • ■ Ja{o)=M)' e M^^^^ and Q' is the transpose of 
Q. 



Remark 3.1 Notice that (Pt{0) = 1 and (pT{—i) = E Q^j=i f^^^^'^i . Furthermore, if fj,{a) = fi, 
then (^5"(— i) = e^-^ since Yl^i ^« — ^- 



Fourier methods for switching jump-diffusions and FSO pricing 



9 



More generally, we get from ([9]) and ()12p 

ift^Tiz) = Et[e''^''T] = (l' ■ e^Q'+' diag(^i(^),...,^M(^)))(T-t) . 

(13) 

= Ejlii«w=,qf (^), qf (-) = Ef=i (e^'^'^' diag(.,(.),...,.„(.)))(T-.))^^ 

for any t € [0,T), being the conditional expectation up to time t. Notice that the charac- 
teristic function of X{T) and Xt^T depends on the state of the Markov chain a(0) and a{t), 
respectively. 

The conditions for applying the transform method both in log-price and log-strike depend on 
the properties of the GTF of X{T), which in turn depend on that of (j)j{z) through Proposition 
13.11 In general, these functions are well defined (and analytic) in some strips of the complex 
plane 

Sj = {zeC: E2[e^(^)T{^(j')j)] < oo}, j = l,...,M. 

Let us define the matrix A{z) = Q' + i diag{'&i(z), . . . ,'dMiz)): clearly the elements of 
A{z)^,n = 1,2,... are polynomials in the ■dj{zys and therefore these are well defined in the 
intersection of the Sj, j = 1, . . . ,M. From the properties of the matrix exponential function 
gA{z) _ ~^^^r~ since the GTF of X{T) is a linear combination of its elements, it im- 

mediately follows that (jl2p and (jl3p are well defined in Hjli '^j ^'^d consequently the transform 
methods can be applied, provided fljli 'Sj 7^ and the payoffs satisfy the proper conditions. 

Remark 3.2 If we set /i(i) = n,cr{i) = a, A(i) = A and (j)i{z) = (j){z) we have that '&i{z) = 0, 
i = 1, ... ,m — 1, and the term exp(Q'T) is the transpose of the transition semi-group of the 
Markov chain. Under these choices we are implicitly assuming a unique regime and eq. U3\} 
becomes the well-known characteristic function of the (single-regime) jump-diffusion dynamic 
(0), (friz) = exp(ze + ^i^V2-iA((/.(z)-l)). This is because I'-e^^'+i diag(^i(^),...,^M-iW,o)){T-t) . 
I{t) = 1' • e*5 C^-*) . i(^t) = X^i^i Ia(t)=j = 1- Hence, with simple linear constraints on the full 
parameter set of our dynamic (B^ we can recover several models: 

1. Black & Scholes model (BS): fii = r, ai = a > 0, Xi = (we consequently set to zero the 
jump variables Y{i) ), i = 1, . . . , M; 

2. Black & Scholes with regime switching model (RSBS): Hi G M, cTj > 0, Qij > 0,i ^ j, 
Xi = (V{i) = 0), i = l,...,M; 

3. Merton jump-diffusion model (JDM): ij,i = r,ai = a>0, Xi = X>0 and the parameters 
of the jump variables Yi = Y , i = 1, . . . , M; 

4- Merton jump- diffusion model with regime switching (RSJDM): fii E M, fij > 0, Qij > 0,i ^ 
j, Aj > and the parameters of the jump variables Y(i) for each regime, i = 1, . . . , M . 

From a computational viewpoint, for a fixed complex z the calculation of (Pt{z) requires the 
following steps: 
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1. calculate ^j{z),j = 1, . . . ,M (eq. (fTTD): 

2. form the matrix A{z) = Q' + i diag(i?i(z), . . . , i?j\,f (-2)); 

3. calculate the matrix exponential ^{z) = exp{A{z) T); 

4. for each starting state of the chain a(0) = j, j = 1, . . . , M, calculate qj{z) = Ylh=i ^kj{z)- 

For M > 2 the cumbersome task is the calculation of $ for which efficient numerical tech- 
niques are available (see Higham (2009)). 

The case M=2. In this case it is possible to give a closed form solution to the matrix expo- 
nential, therefore obtaining an easy-to-implement formula for the characteristic function. The 
following result can be proved either by solving a couple of ODE, as in Buffington and Elliott 
(2002) - Appendix 1, or through a Laplace Transform - based technique, as in Liu et al. (2006). 

Proposition 3.2 Letyi2 he the solutions of the quadratic equation 'ip' + [qi + q2 — '\0)y — \9q2 = 
and 

qi'^W = ^(e^^^^"*H?/i + ^i + 92)-ei'^(^-*)(y2 + gi + g2)) 

qf(^) = ^,{^'''^'^-'\yi + qi + q2--^e)-M^-'\y2 + qi+q2--^e)). 

Then 

= Mt)=i^f{o) + Mt)=2^f{o). 

□ 

It is easy to prove that the functions q^'"^ and ({^ are invariant by changing the order of the 
roots y\ and 7/2- The characteristic function follows from the proof of Prop. 13.11 (see (j32p ): 

^tA^) = e'^^(^)(^-*) (l,(,)=iqf (0(z)) + I,(t)=2qf (^(^))) • (14) 

1 2 l2 

Example 3.1 In our regime switching version of the Merton model we have 4'i{z) = Q^^°-i~2^ i ^ 
i = 1, . . . , M. Then, from ill]) 

§,{z) = zii + hz^aj - iA,(e^^'^»-^^'''? - 1), i = 1, . . . , M. (15) 

It follows that in such a case the characteristic function (Pt{z) is well defined for all z €C In 
the two state model the GFT is easily obtained from 

Some examples of payoff transforms for the typical claims are recalled in Table [H In view of 
our next applications, we show in some details how to get the price of call and put options both 
in log-price and in log-strike transform. Pricing formulas for the other payoffs are reported in 
Table El 
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Call/Put value in log-price transform. From formula (|10p and Table [T] we get for the call 
option 



Co 



P{0,T) 
2tt 



:,fc(i2:+l) 



-dz, U > 1, 



iz — z^ 



(16) 



27r 



log(so)+fc(l-!y) 



+ 00 



'iu(log(sQ)-k) _ 



lPt{—u — iv) 



27r 



+ 00 



-in log(so/-ft')_ 



— — ti^ + iti(l — 
iPt{—u — w) 



-du 



-du, 



3 1/2 — z/ — li^ + m{l — 2v) 

provided the characteristic function evaluated in the integral ()16p is well defined for z S C 
such that '^{z) > 1. By switching from '^{z) > 1 to < we get the value for the put 

option: notice that the put-call parity relation is recovered by moving the integration contour. 
As a matter of fact, alternative formulas can be derived by using residue calculus (see e.g. 
Lewis(2002)), under the proper conditions for ipT{z). The GFT of this payoff has two simple 
poles at z = and z = i with residue — and ^o*^^^ ')^ ^ respectively: by moving the integration 
contour and since the integral must be real, we obtain the following general formula in which 
we stress the dependence on sq, ckq and K: 



Co{so,ao,K) = P{0,T) [R^ + ^ 



iu+oo k(iz+l) ^ 

e-'"'°g(^°VT(-^) _ ^ dz 



(17) 



P{0,T) R^ + -s^qK 
vr 



-\u \og(so/K) _ 



lPt{—u — \v) 



M 



/ . -Q(0)=i 



1/ - n2 + 1-^(1 - 2v) 
q.' {-u-\v) 



du 



-iulog{so/K)0,T 



v"^ — 1/ — -\- 'm{l — 2iy) 

V > 1 

V = 1 
< 1/ < 1 



du 



(18) 





So 

so<^T(-i) 

soV3T(-i) - T ^ = 
so93t(— i) — e'^ V <^ 

where (^t(— i) = E[e^»=i ^'"^'] according to Remark (13. ip and the functions q^'"^(-) are defined in 
{[13D. 



Call/Put value in log-strike transform. As before, if the GFT (/'j(-) are well defined func- 
tions in a properly defined strip of C, from formula (llOp and Tabled! we get for the call option 

from which 

P(0,r) _i^fcei°s(^«)(^+i^)(^T(^-i) 



Co = / e-'^'^- : ^ '-dz u<0. (19) 
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2tt 



Ju{k~log{so)) _ 



(friu + i{v - 1)) 



-iu \og{K/so)) _ 



— V — + in(l 
ifT{u + \{u - 1)) 



2v) 
-du. 



du 



27r ' u'^ -u -V? + m{l-2v) 

The value for the put option and the related put-call parity are obtained again by moving 
the integration contour. Since the residues at the poles z = and z = i of the integrand are 
^oyrC-i) g^j^j -gfc respectively, the application of the residue Theorem gives the following general 
formula for the call price in our RSJD model: 



/ 1 r'w+co 

Co(so, ao, K) = P(0, T)\r, + — I e 



\z — 



-dz 



(20) 



P{0,T)(R,+ 



M 



-iu \og{K/so)) _ 



ipriu + i(i/ - 1)) 



P(0, T)\r, + -sl^'K^ ^°(o)=J- ^ 



v"^ — V — -\- m{l — 2v) 
%''^{u + \{v - 1)) 



-\u\og{K/so)0,T 



v'^ — V — + m{l — 2u) 



du 



du 



(21) 





So 

SQipri-'-i) 



1/ < 
1/ = 
< 1/ < 1 

so(/?T(-i) - ^ = 1 
so</?t(— i) — e'^ V > 1 



Remark 3.3 Let us notice that due to the symmetry of the call payoff, the two approaches give 
in general very similar pricing formulas: in particular from f7^) and / f7P|] it follows that by 
changing z with i — z we can switch from one representation to the other. 



Application of the FFT algorithm. As it is widely known, the transform method deserves 
for an efficient evaluation of derivative prices by means of the FFT algorithm for a proper range 
of the trigger parameter. Actually, if only one option price has to be evaluated for a fixed k, 
there is no need to use FFT. This technique involves two steps: 

1. a numerical quadrature scheme to approximate the integral appearing in the pricing for- 
mula, that we write as 

m = - 

through a A'^-point sum. By using an equispaced grid {un}n=i,...,N of the line {z = u + iv 
C : u G M"^, u = ly} with spacing A, we have 



r+oo 

/ 3? 


e-™'=F(n)" 


du 


Jo 







I{k) « T.N{k) = 
where Wn are the integration weights; 



A 

vr 



N 



-iu„k 



F{Un)Wn 



n=l 



2. given a properly spaced grid of triggers km = ki -\-^{m — 1), m = 1, . . . N , the sum Ti]y{k) 
is written as a discrete Fourier transform (DFT), so that the FFT algorithm can be used. 
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Payoff 


Option value in log-price transform 




P(0,r)i.o^e'=(i-)E-liW)=,/o'"°^» 


(,-m(log(so)-fc)qO,T^_^_j^^ 

!^-^— i^— n'^+m(l— 2i/) 


(iii, z^ < 




P(0, T)i4e('^-)^/'' E -li I.(0)=, lo^ 3? 


"g-iti{log(so)-fe/6)qO.T(„u„ij,) 
I/— a— m 


du, V > a 


min(e^, e'^) 


p(o, r)i.[;e'=(i-) EfLi Uo)=j /o'"°° 3^ 


"e-i«(log(so)-fc)qO,T^_^_j^^- 


du, < < 1 




Option value in log-strike transform 


(efc _ ey)+ 


P(0,r)i4-^e'=-E^iIa(o)=,/o'°°3? 


"e-'"('=-'°s(=o))qO'^(u+i(;y-l))" 
i/-^— u'^+m(l— 2;/) 


du, u > 1 




P(0, r)i.r'''^e-'= 1.(0)=, lo^^ K 


g-m(fc-blog(so))qO.T'(a_i,b+iub) 


du, 1/ > 


min(e^, e'^) 




du, < < 1 



Table 2: Option values for some typical payoffs under the RSJD model. See the Appendix for 
a sketch of their derivation. 



A numerical example. In order to asses the performances of the pricing formulas we consider 
the basic models in Remark (j3.2p . Example (j2.ip . in which the regime switching behavior is 
driven by a two-state Markov chain. We fit these models on a set of observed call prices on the 
S&P 500 index as quoted on March 31, 2009 to get realistic values for the parameters. In the 
data set used for calibrating the models there are 128 call option prices with maturities and strike 
prices ranging from 31 to 272 days and from 525 to 1200, respectively. The value of the index is 
So = 753.89 and the moneyness sq/K ranges from 0.6282 to 1.4360. The average of the bid and 
ask Treasury bill discounts, as available from the Wall Street Journal, were used and converted 
to annualized risk-free rates. The dividend rate q was estimated from the data: in particular 
we used a non linear least squares algorithm which minimize the difference between observed 
call prices and the corresponding Black & Scholes prices evaluated through the available implied 
volatility, constrained to satisfy the put-call parity relations. This procedure was repeated for 
each maturity giving a mean value q = 0.0157 with standard deviation 0.003. 

The numerical implementation was developed in the MatLab® environment. Quadrature 
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, ' ^0.1508 
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20 



25 



30 



Figure 2: Real part of the integrand for a call option for different values of the moneyncss. The plots 
related to our data (0.6282 and 1.4360) are very similar. 



algorithms are needed to evaluate the option prices from (jlSp : adaptive Simpson and Gauss- 
Lobatto quadrature rules, as available in MatLab, performed equally well, for typical values 
of the parameters. As a matter of fact the integrands are not rapidly oscillating and decrease 
sufficiently fast (e.g. see Fig. ([2])). The FFT algorithm was implemented following Lee (2004), 
i.e. by sampling F at the midpoints of intervals of length A, Un = {n — ^)A, n = 1, . . . , and 
taking 7A = 2t:/N. We get 



S7v(fem) = 

vr 



■ TV 
.n=l 



-i(n-i)A(fci+7(m-l))^(-(,^ 



1 



A) 



vr 



N 

Ve"'^^("" 



l)(m- 



n=l 



where f{n) = F{{n — i)A)e 



-i(n-i)Afci 



In this case we used 7 = 0.01 and ly = 0.5. 



For the calibration we minimized the sum of squared errors by using the constrained min- 
imization routine in MatLab. In fact, for the regime switching models we have to add the 
constraint ai > o"2. The results obtained are reported in Table ([3]): RMSE and relative errors 
^~^o^o,K) calculated in the four cases. In Table (HD we report out-of-sample performances 
of each fitted model: these were obtained by calculating the deviation from five call option prices 
having a much longer maturity, i.e. 631 days and moneyness ranging from 0.7539 to 1.0052. 
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Figure 3: Implied volatility calibration. 



4 On the pricing of forward starting options 

Forward starting options are well-known exotic derivatives, depending on an underlying asset 
characterized by the payoff 

Ut{S{T), k) = S{T) - KS{t*), (22) 

where t* € (0, T) is the determination time and k € (0, 1) is a given percentage. They are the 
building blocks of the so-called cliquet options and are used in many different context. 

In this Section we provide a simple valuation formula for the price at time i = of this claim 
where the underlying S{t) follows the regime-switching jump diffusion dynamic introduced in 
Sect. 2. Furthermore, we assume that fi{a) = r, the risk-free rate, in such a way V = Q. The 
risk-neutral price is therefore given by 

Uoiso, ao, k) = E[e-'-^(5(r) - ))+] = E[e-^'^(soe^(^) - KS{t*))+]. 

Notice that in general from the determination time t* on, the price is equal to that of a 
standard call option, being the strike a known constant. By denoting with Ei[-] the conditional 
expectation w.r.t. information up to time t, J-f, we have 

Ct{S{t),a{t),K) = Et[e-'^^-'HS{T)-K)+] = e-^(^-*)Ei[(5(t)e^*.^ - 
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BS 


RSBS 


JDM 


RSJDM 


0-1 


0.3645 


0.4462 


0.1341 


0.2725 


0-2 




0.3296 




0.1350 


Ai 






7.9958 


6.8393 


As 








0.8590 


ai 






-0.1280 


-0.1398 


0.2 








-0.3423 


h 






0.0011 


0.0877 










0.1593 






9.6199 




6.5075 


92 




0.0002 




0.0020 


RMSE 


4.6947 


3.9177 (io = 1) 


3.8715 


0.6126 {io = 1) 


Rel. err. 


(-0.0242,0.1398) 


(-0.0141, 0.0631) 


(-0.0508, 0.1547) 


(-0.0080, 0.0191) 


0.0888 


0.0353 


0.0675 


0.0041 



Table 3: Implied parameters and in-sample calibration performances of the models. The mon- 
eyness of these options ranges from 0.6282 to 1.4360. In the last rows we report the range and 
the mean of the relative pricing error ^"^^^ ■ 





BS 


RSBS 


JDM 


RSJDM 


RMSE 


14.8555 


6.4732 


17.1020 


5.1116 


Mean Rel Err. 


-0.1852 


-0.0631 


-0.2142 


-0.0634 



Table 4: Out-of-sample performance of the models. The moneyness of these options ranges from 
0.7539 to 1.0052. 

Therefore, \i K = nS{t) we get that 

Ct{S{t),a{t), KS{t)) = S{t)Ct{l, a{t),K). 

Hence, by the law of iterated conditional expectations, 

no(so,ao,'^) = E[e-''^(5(t*)e^'*.^ - tiS{t*))+] = Ct*{S{t*),a{t*),KS{t*))] = 

= E[5(r)e-'^**a*(l,a(r),Ac)]. 
Since S{t)e~^^ = S{t)/B{t) is a Q-martingale we can introduce an equivalent measure as 

Hi) = ^«!|, = E[L(T)1 = 1, 

^ ' dQ ^ B{t) so ^ ^ 

from which we get 

E[5(r)e-'-**a.(l,a(r),K)] = soE«" [a* (l, «(**), /^)]. 

Notice that this property is fairly general: in fact a{t) is not restricted to be a Markov chain. 
On the other hand, in our model we don't need to further specify the Q'^ dynamic of the price 
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process, since the Markov chain is not affected by this change of measure. Hence, since the chain 
is assumed to be stationary, by denoting with vr = (tti, . . . , ttm)' its invariant probabihty, we get 

M 

Uo{so,k) = soE'^'[Ct*{l,a{t*),K)] = sq vr,a* (1, J». 

j=i 

The price of the forward starting option is therefore the mixture of cah option prices evaluated 
at the determination time t* under each regime weighted by the corresponding probabihty. 
By using transform representation for the cah option value, e.g. (jlSp . we get the following 
proposition. 

Proposition 4.1 Let the underlying be characterized by the SDE Then the price at time 
t = of the Forward Starting Option Ii22\) with determination time t* and percentage k € (0, 1) 
is given by 



M / i_„ ,4-^ r Q-iu\og{K) t* ,T , _ ■^ 



R, + / 9? 



■J 



u"^ — m{l — 2iy) 



du 



where {TTj}j=i^,„^m is the stationary probability of the Markov chain, v and Ry following from 



As a byproduct of the last proposition, by restricting our model to a unique regime (see 
Remark (j3.2p ). we get a simple formula for pricing a FSO in a Levy model with finite activity. 

The impact of model choice on the prices of the Forward Starting options is shown in figures 
dH) , ([5]) and ([6]) as a function of the determination time for three different values of the percentage 
K. The parameters of each model are those estimated in our numerical example (Table ([3])). 



Pricing FSO in a general stochastic volatility model. In Chourdakis (2004) a regime- 
switching diffusion was considered to approximate a general stochastic volatility model 

dX{t) = fi{v{t))dt + a{v{t))dW{t) + [ yp'^{dy,dt), (23) 
dv{t) = a{v{t))dt + b{v{t))dZ{t), (24) 

where X{t) = log(5(t)) and p°'{dy,dt) has intensity X{t,v,dy) = X{v)my{y)dy, m„ being the 
probability measure which characterizes the jump component. Then, under some conditions 
on the coefficients a(-) and 6(-), the diffusion process ()24p can be approximated by a finite 
state Markov chain defined on a grid which is the discretization of the domain of v. The 
approximating scheme defines a generator Q*^ for the Markov chain depending on e and on the 
functions a(-) and b{-) evaluated at the points of the grid = {wf, . . . ,vl^j}. As reported in 
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Chourdakis (2004) the generator Q"^ = {qfj}ij=i,...,M where 

^bHv'^) - ^X^'^), i = j-l 

^b'^ivj), i = j 

(25) 

^b\v^^) + ^^a{v^^), i = j + l 

produces accurate results for coarse volatihty grids. The resulting approximated process X^{t) 
follows therefore a RSJD dynamic and its characteristic function is obtained from Prop. (j3.ip . 
Correspondingly, option prices can be calculated by means of the Fourier transform techniques 
presented in Sect. 3. Convergence properties as well as computational considerations as e — > 
are discussed in Chourdakis (2004) where a number of cases are studied. 

This technique combined with our Proposition (j4.ip suggests the following scheme for pricing 
FSO under a general SV model: 

1. approximate the model with a regime-switching diffusion X''{t): this amounts to build the 
generator of the Markov chain, defined by ([25]); 

2. evaluate call option prices C| at time t* under each regime j = 1, . . . , M through formula 
([H]) or din) with S{t*) = 1 and K = k; 

3. calculate the price no(so, k) = sq Z^jli '^j^j where the coefficients VTj's are the solution of 
ttQ' = 0. 
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5 Conclusion 

In this paper we considered the problem of valuing the price of a European contingent claim 
when the underlying dynamic follows a Levy process of I type whose parameters are modulated 
by a continuous time and finite state Markov chain. These kind of processes are known to cap- 
ture specific features of financial time series, such as volatility clustering and structural breaks. 
On the other hand, they can equally be used to approximate very general stochastic volatility 
processes. Following the well established relationship between option prices and Fourier trans- 
forms, we obtained almost closed- form solutions (up to a numerical integration) for European 
style options, both in log-price and in log-strike space. An example of calibration for the regime- 
switching version of the Merton jump-diffusion model is also presented for a daily set of call 
option data on the S&P 500. 

Furthermore, as a practical application of the Fourier transform methodology we obtained an 
almost closed-form solution to the problem of valuing a Forward Starting option in our general 
regime-switching jump-diffusion dynamic. This result can be jointly used with the approximation 
scheme of stochastic volatility models to get a feasible algorithm for FSO pricing under a very 
general dynamic. 

6 Appendix 

Proof of 12.11 Let us define the occupation times for the Markov chain a{t) in [0,T], Tj = 
/q ^a{s)=ids, i = 1, . . . , M. We immediately have that ^ T, = T. Now, given a sample path 
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Figure 6: Forward starting option prices. 



of the chain a{t), < t < T, we can define 

Ai = (J [T,,r,+i) (26) 

£: 0!(t)=i,Ti<t<Ti^i 

^^(A^) = Yl {Nn+^-Nr,) (27) 

£: a{t)=i,Ti<t<Ti+i 

Z(A,) = {W{n+i)-W{n)). (28) 

t: a{t)=i,Ti<t<Ti+i 

Since each Aj is the union of non overlapping intervals, the corresponding random variables 
A^(Aj) and 2'(Aj) are distributed as a Poisson variable Poiss(Airj) and as a Normal variable 
M{Q,Ti), respectively. Furthermore, iV(Aj) ± iV(Aj) and Z(Ai) _L Z(Ai), for i / Jl. By 

denoting with the /c-th jump magnitude relative to regime i, we have 

j=l fc=l 1=1 i=l 

By defining ^(a) = fj,{a) — ^cr^(a) — X{a)K{a) and 

Hr(ri, . . . , T^) = r ^(a(s))ds = ^(i)r, 



■^Here, X J-Y means that X and Y are independent. 
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then X{T) admits the following representation: 



M N{Ai) 



X{T) = ~T{Ti,...,T^) + Y,<^{i)Z{^i) + Y. E 



(30) 



i=l k=l 



Proof of 13.11 Let (l)j{z) = 'E[e^^'^^'^^^^'^^] be the generalized Fourier transform of the jump 
magnitude. From the representation ()30p we can easily calculate the characteristic function of 
X{T), conditional to J^^: 

E[ei.x(T)|^aj ^gi.=T(Ti,...,T„) 5,jgi.Er=i-0-)^(A.)|j-a] E[e'^ ^^'=1 S^=^'' '^'^ | J?] . 

The first expected value is simply obtained as 

while the second, since A^(Aj) ± N{Aj) for i ^ j, is 



Finally, we have 

ifriz) = E 



(31) 



Actually, the exponent in (|31|) is a linear function of the sojourn times Ti, • • • ? -^m? 

the 

characteristic function of which are well-known. As a matter of fact, we have 



ipriz) = E 



e ^3 



(32) 



where ??i(2;) = ^^^(a:) — 'i?m(-z)) being = T — (Ti + . . . + T^-i). Since it can be proved (see e.g. 
Buffington and Elhott (2002)), that 



E 



e 



-,_/ . gQ'+idiag(^i,...,,?„_i,0)T . j(^Q^^ 



formula (jl2p follows, the second equality being a consequence of the property of matrix expo- 
nential exp(0) exp(A) = exp(^/ + A). 



Derivation of Tables [2l Payoff e'^^I^x^K- From the third row of Tableland formula (fTO 
we get for log-price transform 



no/P(o,r) 



1 



{a+iz)k/b 



2vr Jiu-oo CL + iz 



-ipx{—z)dz 



2tt J-oo V — a — m 
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~iu{\og{so)-k/b) 



ifxi—u — iv) 



V — a — m 



du, V > a; 



for log-strike transform 



^\^X{T)+\og{so){k)] = E[- 



,{a+izb)(\og{so)+X(T)) 



^{a+izb) log(so) 



IZ 



IZ 



-ipria + izb) 



from which 



no/p(o,r) = — 



Ja+izk) log(so) 
-izk "= 



27r 



a—bui^k 



vr 



+00 



11^— 00 

+00 g-m(fc-61og(so)) 

-00 in - 1/ 
g-m(fc-61og(so)) 



ipT{cL + \hz)dz 



-LpT{a — vb + ibu)du 



m — V 



iprio- — ub + ibu) 



du, 1^ > 0. 



Payoff min(e^,e'^). From the fourth row of Table [T] and formula (|1U|) we get for log-price 
transform 



no/p(o,r) 



1 

2^ 



ii^+oo 



-iz log(so) 



,fc(i^+l) 



Z^ — iz 



ipT{—z)dz 



27r 



+00 



-m(log{so)-fc) 



vr 



SqC 



+00 



^2 _ ^2 _j_ _j_ [y^(2iy — 1) 
g-i«(log(so)-fc) 



ipT{—u — \v)du 



V? — v"^ + V + iii(2i/ — 1) 



(Pt{—u — iu) 



du, < ly < 1; 



for log-strike transform 



Efn 



X(T)+log(so 



,{l+iz){log(so)+X(T)) 



^{1+iz) log(so) 



IZ 



Z^ — iz 



-lPt{z - i) 



from which 



no/P(o,T) 



2tt 



iu+oo 



. , e(l+iz)log(so) 
-\zk „—izk ^ 



— IZ 



(Pt{z — i)dz 



vk „l—u 



27r 



-e s 







-iM(fc-log(so)) 



vr 



e 6q 



5? 



^2 _ j^2 _j_ _j_ i^(2z/ — 1) 

-m(fc-log(so)) 



(/5T(n -|- i(i^ — l))du 



^2 _ j^2 _|_ _j_ 'iii{2v — 1) 



(/9r(n + - 1)) 



du, < < 1. 



By substituting in the previous formulas the GFT (|13p of the RSJD model, we immediately 
get the entries of Table [21 



Fourier methods for switching jump-diffusions and FSO pricing 



23 



Acknowledgments 

The financial support of tlie Research Grant: PRIN 2008, Probability and Finance, Prot. 
2OO8YYYBE4, is gratefully acknowledged. Moreover, the author would like to thank the partic- 
ipants to the workshop " Stochastic Volatility, Affine Models and Transform Methods" organized 
by Prof. S. Herzel at the School of Economics of the University of Roma - Tor Vergata, 15-16 
April 2010. 



References 

Bakshi G. and Madan D. (2000), Spanning and derivative-securities valuation. Journal of 
Financial Economics, 55, pp. 205-238. 

F. Biagini, Y. Bregman, and T.Meyer-Brandis, Pricing of catastrophe insurance options written 
on a loss index with reestimation. Insurance Math. Econom., 43 (2008), pp. 214-222. 

Bollen N. P. B. (1998), Valuing options in regime-switching models. Journal of Derivatives, 6, 
pp. 38-49. 

K. Borovkov and A. Novikov, On a new approach to calculating expectations for option pricing, 
J. Appl. Probab., 39 (2002), pp. 889-895. 

Boyarchenko S., Levendorskii S., American options in regime-switching models. SIAM J. Con- 
trol Optim. 48 (2009), no. 3, pp. 1353-1376. 

Boyle P. and Draviam T. (2007), Pricing exotic options under regime switching. Insurance: 
Mathematics and Economics, 40 , pp. 267-282. 

Buffington J. and Elliott R. J. (2002), American options with regime switching. International 
Journal of Theoretical and Applied Finance, 5 , pp. 497-514. 

Carr P. and Madan D.B. (1999), Option valuation using the Fast Fourier Transform, Journal 
of Computational Finance, 2, pp. 61-73. 

Cherubini U., Delia Lunga G., Mulinacci S., Rossi P., Fourier Transform Methods in Finance, 
Wiley, 2009. 

Chourdakis K. (2004), Non- Affine option pricing. The Journal of derivatives, 2004, pp. 10-25. 

Chourdakis K. (2007), Levy process driven by stochastic volatility, Asia-Pacific Financial Mar- 
kets, Vol. 12, No. 4, pp. 333-352. 

Cont R. and Tankov P. (2004), Financial Modelhng with Jump Processes, Chapman & Hall/CRC 
Financial Mathematics Series. 

Di Graziano G. and Rogers L. C. G. (2009), Equity with Markov-modulated dividends, Quan- 
titative Finance, Vol. 9, No. 1, pp. 19-26. 



Fourier methods for switching jump-diffusions and FSO pricing 



24 



Di Masi, G.B., Kabanov, Y.M., Runggaldier,W.J. (1994), Mean-variance hedging of options on 
stocks with Markov volatihty, Theory of Probability and Its Applications, 39, pp. 173-181. 

Duan, J.C., Popova, I., Ritchken, P. (2002), Option pricing under regime switching. Quantita- 
tive Finance, 2, pp. 1-17. 

D. Dufresne, J. Garrido, and M. Morales, Fourier inversion formulas in option pricing and 

insurance, Methodol. Comput. Appl. Probab., 11 (2009), pp. 359-383. 

E. Eberlein, K. Glau, and A. Papapantoleon, Analysis of Fourier transform valuation formulas 

and applications, Appl. Math. Finance, 17 (2010), pp. 211-240. 

Edwards C. (2005), Derivative Pricing Models with Regime Switching. A General Approach, 
The Journal of Derivatives, Vol. 13, No. 1, pp. 41-47. 

Elliott R. J. and Osakwe G. J. U. (2006), Option pricing for pure jump processes with Markov 
switching compensators. Finance and Stochastics, 10, pp. 250-275. 

Elliott R. J. , Siu T. K. and Chan L. (2007), Pricing volatility swaps under Heston's stochastic 
volatility model with regime switching. Applied Mathematical Finance, Vol 14, 1, pp. 41- 
62. 

Guo X., (2001), Information and option pricing Quantitative Finance, 1, pp. 38-44. 

Guo X. and Zhang Q. Z (2004), Glosed-form solutions for perpetual American put options with 
regime switching, SIAM Journal of Applied Mathematics, 64 , pp. 2034-2049. 

Hamilton, J.D., (1989), A new approach to the economic analysis of non-stationary time series, 
Econometrica, 57, pp. 357-384. 

Hamilton, J.D., (1990), Analysis of time series subject to changes in regime. Journal of Econo- 
metrics, 45, pp. 39-70. 

Hardy, M.R., (2001), A regime switching model of a long term stock-returns. North American 
Actuarial Journal, 3, pp. 185-211. 

Heston, S. L. A. (1993), Glosed-Form Solution for Options with Stochastic Volatility with 
Applications to Bond and Currency Options, Rev. Fin. Studies, vol. 6, 327-343. 

Higham N. J., (2009), The scaling and Squaring Method for the Matrix Exponential Revisited, 
SIAM Review, Vol. 51, No. 4, pp. 747-764. 

F. Hubalek, J. Kallsen, and L. Krawczyk, Variance-optimal hedging for processes with station- 

ary independent increments, Ann. Appl. Probab., 16 (2006), pp. 853-885. 

T. R. Hurd and Z. Zhou, A Fourier transform method for spread option pricing, SIAM J. 
Financial Math., 1 (2010), pp. 142-157. 

Jiang Z., Pistorius M. R. On perpetual American put valuation and first-passage in a regime- 
switching model with jumps. Finance Stoch. 12 (2008), no. 3, 331-355. 



Fourier methods for switching jump-diffusions and FSO pricing 



25 



Jobert A. and Rogers L. C. G. (2006), Option pricing with Markov- Modulated dynamics, SIAM 
Journal of Control and Optimization, Vol. 44, No. 6, pp. 2063-2078. 

Khaliq, A. Q. M.; Liu, R. H. New numerical scheme for pricing American option with regime- 
switching. Int. J. Theor. Appl. Finance 12 (2009), no. 3, 319-340. 

Konikov M. and Madan D. B. (2002), Option Pricing Using Variance Gamma Markov Chains, 
Review of Derivatives Research, Vol. 5, No. 1, pp. 81-115 

Kruse S. and Nogel U. (2005), On the pricing of forward starting options in Heston's model on 
stochastic volatility. Finance and Stochastics, 9, pp. 233-250. 

Lee R. W. (2004), Option pricing by transform methods: extensions, unifications and error 
control, Journal of Computational Finance, 7, pp. 51-86. 

Lewis A.L. (2002), A simple option formula for general jump-diffusion and other exponential 
Levy processes, Working Paper, Optioncitynet.net. 

Liu, R. H. Regime-switching recombining tree for option pricing. Int. J. Theor. Appl. Finance 
13 (2010), no. 3, 479-499. 

Liu R. H. , Zhang Q. , and Yin G. (2006), Option pricing in a regime-switching model using 
the fast Fourier transform. J. Appl. Math. Stoch. Anal., Vol. 22, pp. 1-22. 

R. Lord (2008), Efficient pricing algorithms for exotic derivatives, PhD thesis, Univ. Rotter- 
dam. 

Naik, V. (1993), Option valuation and hedging strategies with jumps in the volatility of asset 
returns. Journal of Finance, 48, pp. 1969-1984. 

S. Raible (2000), Levy processes in finance: Theory, numerics, and empirical facts., PhD thesis, 
Univ. Freiburg. 

Ramponi A. (2009), Mixture Dynamics and Regime Switching Diffusions with Application to 
Option Pricing. Methodol. Comput. Appl. Probab., DOI 10.1007/sll009-009-9155-l 

Runggaldier W.J. (2003) . Jump-Diffusion models, in "Handbook of Heavy Tailed Distributions 
in Finance" (S.T. Rachev, ed.). Handbooks in Finance, Book 1 (W.Ziemba Series Ed.), 
Elesevier /North-Holland 2003, pp. 169-209. 

Timmermann A. (2000), Moments of Markov switching models. Journal of Econometrics, 96, 
pp. 75-111. 

Yao D. D., Zhang Q. and Zhou X. Y. (2006), A regime-switching model for European options. 
Stochastic Processes, Optimization, and Control Theory Applications in Financial Engi- 
neering, Queueing Networks, and Manufacturing Systems (H. M. Yan, G. Yin, and Q. 
Zhang, eds.). Springer, New York, pp. 281-300. 



